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Abstract 

Fermi 7-ray telescope data of GRB 080916C with ^ 10^^ erg in apparent isotropic 
7-ray energy, show a several second delay between the rise of 100 MeV - GeV radiation 
compared with keV - MeV radiation. Here we show that synchrotron radiation from 
cosmic ray protons accelerated in GRBs, delayed by the proton synchrotron cooling 
timescale in a jet of magnetically-dominated shocked plasma moving at highly relativis- 
tic speeds with bulk Lorentz factor F ~ 500, could explain this result. A second gen- 
eration electron synchrotron component from attenuated proton synchrotron radiation 
makes enhanced soft X-ray to MeV 7-ray emission. Long GRBs with narrow, energetic 
jets accelerating particles to ultra-high energies could explain the Auger observations of 
UHE cosmic rays from sources within 100 Mpc for nano-Gauss intergalactic magnetic 
fields. The total energy requirements in a proton synchrotron model are oc F^^/'^. This 
model for GRB 080916C is only plausible if F < 500 and the jet opening angle is ^ 1°. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mechanisms: 
nonthermal 

1 Introduction 

An integrated fluence of 2.4 x 10~^ erg cm~^ was measured from GRB 080916C with the 
Large Area Telescope (LAT) and Gamma ray Burst Monitor (GBM) on the Fermi Gamma 
ray Space Telescope, with one third of this energy in the LAT At a redshift z = 
4.35 ± 0.15 [2], GRB 080916C has the largest apparent energy release yet observed from a 
GRB. A significant ~ 4.5 s delay between the onset of > 100 MeV compared to the ~ 8 keV 
- 5 MeV radiation is found (the characteristic duration of the GBM emission is ~ 50 s). 

The spectrum of GRB 0809 16C was fit by the smoothly connected double power-law 
Band function [3] to multi-GeV energies, though with changing Band spectral parameters 
and peak photon energy in different time intervals. The emergence of delayed spectral 
hardening is represented by a Band beta spectral index changing from /3 = —2.6 in the first 
3.6 seconds following the GRB trigger to /? = —2.2 at later times jl]. Here we show that a 
hard spectral component arising from cosmic-ray proton synchrotron radiation explains the 
delayed onset of the LAT emission. If GRBs accelerate UHECRs, then the delayed onset 
of the LAT emission after the GBM trigger should be a regular feature of GRB spectral 
evolution. 
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2 Proton acceleration and radiation 



GRB blast wave calculations usually treat electrons Ul |^ , but protons and ions will also be 
accelerated if they are present in the relativistic flows in black-hole jet systems. Here we 
consider protons accelerated in GRB blast waves to such energies that they can efficiently 
radiate hard ~ GeV - TeV photons by the proton synchrotron mechanism [HI EJ El El [IHl 
lllj . The highest energy photons are reprocessed by 77 — )• e~^e~ opacity |12j to make an 
injection source of electrons and positrons that cool by emitting < GeV electron synchrotron 
radiation. Two delays arise, the first from the time it takes to accelerate protons to a 
saturation Lorentz factor where the acceleration rate equals the synchrotron rate. A second 
delay arises during which sufficient time passes to build up the spectrum of the primary 
protons so that they become radiatively efficient in the LAT band. Observations of GRB 
080916C are consistent with the delayed onset between GBM and LAT emission being 
caused by the second delay where the evolving proton cooling synchrotron spectrum sweeps 
from higher energies into the LAT waveband. This radiation is emitted from a jet of 
magnetically-dominated shocked plasma with > 1, where is the ratio of magnetic- 
field to proton/particle energy density in the plasma. 

The rapid variability, large apparent luminosity, and detection of high-energy photons 
from GRBs can be understood if this radiation is emitted from jetted plasma moving with 
bulk Lorentz factor L ^ 1 towards us. Detection of 3 GeV and 13 GeV photons from GRB 
080916C suggests T3 = L/IOOO ~ 1 (Ref. [I] and below. Section 3). For variability times 
~ 1 s and instantaneous energy fluxes $ = 10~^$_5 erg cm~^s"^, the internal radiation 
energy density in the fluid is u'^ iircij^^/ {inR'^cT'^) ^ {1 + z)'^d'j^^/ (T^c^tl), implying a 
characteristic jet magnetic field of 



Lgt^ (S 

where primes refer to the comoving frame, pb is the baryon-loading parameter giving the 
relative energy in nonthermal protons compared to 7-rays. The last relation in equation ([T]) 
assumes T rmin from the opacity condition r-yy = 1, which can be written as 

^ r ^Trfi(l + ^)Vgei li/6 

where is the uF^, flux at photon energy e = 2r^/[(l -|- z)^ei] (in mgC^ units) and mgC^ei 
is the highest energy photon. Thus the total jet energy is oc (sPb^j- A characteristic field 

~ 10 - 100 kG is consistent with the absence of a distinct self-Compton synchrotron 
component in GRB 0809 16C which happens when the magnetic energy density is larger 
than the nonthermal electron energy density, or ^^Pb 1- 

In Fermi acceleration scenarios, protons gain energy on timescales exceeding the Lar- 
mor timescale, implying an acceleration rate 7acc,p = (l)~^eB' /rUpC, where i;^"^ ^ 1 is the 
acceleration efficiency. Equating the acceleration rate with the synchrotron loss rate gives 
a saturation Lorentz factor for protons, namely 

, ^mp/B^Y^^ 2 X 108 

where a/ = 1/137 is the fine structure constant, Bcr = rri^c' /eh = 4.4x 10^^ G is the critical 
magnetic field, and the mean magnetic field of the radiating region is B' = lO^iJc G. The 
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observer measures a time 



''sat 



l + zmlc I 67r(j) _ O.OIa/0/10 ^ 



me V earB''^ T^B'^^^ 



for protons to reach ^'sat,p- 

The proton synchrotron saturation frequency (in m(.(? units), corresponding to the 
proton synchrotron frequency of protons with 7^ = ^[at^p^ 

r , T/d) m„ 27 7 Ts 

1 + z ''^ I + z nieoaf 0/10 

This frequency corresponds to a photon energy of ~ 8 TeV for GRB 080916C. Analogous 
to electron blast wave physics [L3\, proton synchrotron losses make a cooling break at 
the proton cooling Lorentz factor 7^^ = "l'sat,p{tsat/t)-, which is obtained by equating the 
synchrotron energy-loss timescale with the comoving time t' = Tt/ {\ + z). Consequently the 
proton synchrotron cooling frequency for protons with 7^ k, 7^^ is £c,p{t) = ^sat.pitsat/t)'^ 
for our simplified model of continuous acceleration and uniform injection of particles during 
the first ~ 8 s to the observer. 

Photons with energies above ~ 1-10 GeV are strongly attenuated through 77 — )• e^e~ 
processes in the source, inducing a nonthermal injection that makes a second-generation 
electron synchrotron component. Proton synchrotron photons radiated by protons with 
I'p — l'sat,p materialize into electrons and positrons with Lorentz factor 7^ w The 
electron synchrotron saturation frequency from this second generation of synchrotron emis- 
sion is observed at 

3 r B' fmp 27 V lO^T^B', 

esat,e - 2 1 + ^'^ B,r\me I6af J ^ (0/10)2 ■ 

Therefore the second-generation synchrotron radiation cuts off above ~ 600 MeV for GRB 
0809 16C. This will be below the LAT sensitivity because only the proton synchrotron ra- 
diation from protons with 7p « 7sai,p contributes to the second-generation electrons and 
positrons that make emission at e ~ £sat,e and the number of such protons is small for a 
steep injection proton spectrum. 

The direct proton synchrotron radiation from cooling protons equals the electron syn- 
chrotron saturation frequency when £c,p{t) = £sat,p{tsat/t)^ = £sat,e, which is observed to 
take place at time 



_ / -Ber- rUe Q4:af _ 4 1 + z rUpcBcT [ m^ , . 

d - sat^(t) ^, mp 81 ~ 3 r e5'2 \j rUe ' ^ ' 

For GRB 080916C, t^i ~ 1.4(0/lO)/r3 ^i^^^^ s, which corresponds (depending precisely on 
£sat,e) to the time required for proton synchrotron radiation to become strong in the LAT 
energy band. We propose this effect as the reason for the delayed onset of the LAT emission. 

In addition to the direct proton synchrotron radiation, pair synchrotron radiation is 
formed from the internal attenuation of proton synchrotron photons to make ultra-relativistic 
pairs. If protons are accelerated with a number index k (i.e., iV(7) oc 7~^), then the cooling 
spectrum breaks from a proton distribution with index k to a steeper one with index A; + 1. 
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The proton synchrotron uF^, flux, in the absence of 77 opacity, is simply 

where fefc^^ = ^/Ob^'cBM, I'minp ~ -'^r'eZ; the relative Lorentz factor of the relativistic wind 
in the jet and the shell material, and 

\ 2 y {£min,p/ £sat,p)^'^ ^^^'^ ~ 1 

* = [ln(7Lt,p/7Ln,p)]"' when A; = 2. 

The direct proton synchrotron flux from GRB 080916C is attenuated at all energies from 
e^^ 1 - 10 GeV to e w £sat,p, equation p), to make a second-generation electron injection 
spectrum with the same form as equation (81), though with subscript p — )• e spectral indices 
{3-k)/2 {3-k)/4 and {2-k)/2 {2-k)/4 and f!:X = \fefaZ- The second-generation 
spectrum has a low-energy cut-off related to e^^, below which it receives no further injection 
pairs. Note that ec,e(i) = ^satA'tsat/t)'^, where esat,e is given by equation 

Fig. 1 shows the evolving proton synchrotron spectrum for constant injection with time 
using the above relations for parameters of GRB 080916C, with the Band function in interval 
(a), the first 3.6 s following the trigger. The Band function in the ~ 10 keV to < 100 MeV 
range can be due to synchrotron radiation by primary electrons, or thermal radiation from a 
jet photosphere. For the proton spectrum, we use k = 2.3, pb = 17, (b = 5, and 7mm p ~ 
In this simple picture of the GRB jet, the protons are accelerated to the saturation Lorentz 
factor within ~ 0.01 s, making a prompt second-generation electron synchrotron spectrum 
too weak to be detected, followed after a few seconds by strong direct proton synchrotron 
emission that sweeps into the LAT band from high energies. The time delay in the high- 
energy 7-ray flux occurs even for variable injection, provided it occurs over ~ 8 s to the 
observer. Although the data for GRB 080916C are consistent with a Band function, spectral 
analyses with the addition of a proton synchrotron component will determine whether this 
more complicated model is compatible with Fermi LAT GRB data. 

Implicit in the treatment is that the wind impacts either a shell ejected earlier by 
the central engine or some roughly uniform density material that already existed in the 
surrounding medium. The extent of this material is (r) ~ cr'^At/{l+z) ~ 5x 10^^r|(At/8 s) 
cm. The magnetic field is thus assumed to be roughly constant over the first ~ 10 s because 
of the constant density medium that is being swept up at the shock. This assumption is 
supported by the roughly constant value of electron synchrotron energy flux $gbm measured 
with the GBM, in the same time interval [1]. 



3 Total Energy of GRB 080916C 

The comoving synchrotron cooling timescale of an ion with atomic mass A and charge Z is 
given by = (74^/Z^)(mp/me)^(67rmec/fTT-B^7). The comoving peak synchrotron photon 
energy is e'^y^ = {Z/A){me/'mp)Bj'^ / Bcr- Equating the observer time tsyn = (1 + z)t'gy^/T 
for the emission to be radiated at measured energy nieC^esyn = ^('syn/i^ + z) = lOO-Eioo 
MeV implies a comoving magnetic field B'{G) = 2.0 x 10^ A^/^Z-''/^E^Q^^Uyri^{s) and an 
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E [GeV] 

Figure 1: Synchrotron model for the 7-ray spectrum of GRB 080916C. The I'Fi, flux in the 
GBM range is produced by primary nonthermal electron synchrotron radiation. The heavy 
dark curves show the fitted Band spectrum for intervals (a), from to 3.6 s after the trigger, 
and interval (b), from 3.6 s to 7.7 s. A strong proton synchrotron component is formed at 
high energies and sweeps into the LAT band after several seconds, thus making the time 
delay between the GBM and LAT emissions. A weak pair synchrotron 7-ray component 
is formed in the GBM and LAT band from second-generation pair synchrotron radiation 
formed by proton synchrotron photons that are attenuated by 77 processes in the jet. Here 
we assume uniform injection over 8 s, and take (p = 10. At a fixed photon energy below 
the cooling frequency, the i^Fiy flux increases oc t for both the proton synchrotron and pair 
synchrotron radiation. 
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isotropic jet power, dominated by magnetic-field energy, given by 



letting the blast-wave radius R ^ r'^ct/{l + z). For Fe {A = 56, Z = 26), the power 
requirements are reduced by a factor ~ 0.17. Here, however, we consider only proton 
synchrotron radiation. 



Eq. (10) shows that Lb oc F^^/s [hJ. The absolute jet power varies as the square of 
jet opening angle 9j. In Ref. [Ij, only the uncertainty in the redshift was used to provide 
uncertainty in T^rnin using simple opacity arguments. Furthermore, validity of the cospa- 
tial assumption that the soft photons are found in the same region as the hard photons was 
assumed in the calculation of Tmin- We use a likelihood ratio test to calculate the chance 
probability to detect a photon with energy greater than the maximum measured photon 
energy in a given time bin. For a —2.2 photon number spectrum, this test gives values of 
F = (0.90, 0.82, 0.75)Fmm for exponential escape and F = (0.79, 0.58, 0.36)Fmin for slab or 
spherical escape at the {la, 2a, 3a) levels, respectively. 

If Type Ib,c supernovae are the progenitors of long GRBs like GRB 080916C, then 
6j > 0.8° [H]. Taking a conservative lower limit = 0.5 gives the absolute energy 
requirements for GRB 080916C of £: ^ 3 x W^^{T3/0.5y^/^{ej /I deg)"^ {tsyn/8 sf^^ erg 



after integrating eq. (10) over time and multiplying by a two-sided jet beaming factor 
fb = 1.5x 10-^{ej/l deg)2. The rotational energy available in a core collapse supernova 
could be as large as f« 5 x 10^^ erg for a 10 Mq core |[15] . 

Although it appears to be a coincidence that F ~ Tmin in our model, which is required 
because of excessive powers when F > F^m, more complicated geometries might relax 
the bulk Lorentz factor requirement further. If the inner engine makes the prompt MeV 
radiation and residual shell collisions at larger radii make LAT 7-ray photons, then F could 
be as low as ~ 300 [16]. In this case, the absolute energy release could be as low as 
2 X 10^^ (^j /I deg)^ erg. Even though the model in Ref. |16| provides a separate explanation 
for the delayed onset (and predicts that the variability timescale of the > 100 MeV radiation 
is longer than the keV/MeV radiation), lack of cospatiality does not guarantee delayed onset. 
For example, photospheric emission with leptonic emission from internal shells would neither 
be cospatial nor necessarily exhibit a delayed onset. 

The jet break time with apparent isotropic energy release ~ 2 x 10^^ ergs is tbr = 
0.3{9j/l deg)^^/^nQ ^^'^ d. The jet break would have taken place before Swift slewed, at 
~ To + 17.0 hr, to GRB 080916C pT], with a hard electron spectrum to explain the shallow 
X-ray decay observed > 10^ s after the GRB [18j. 

A ~ 3 s and ~ 2 s delayed onset of > 100 MeV emission detected in another long GRB 
090902B [19] and in the short GRB 081024B [20], respectively, could be explained with 
proton synchrotron radiation as we discussed here, although lack of redshift measurement 
for GRB 081024B prohibit calculation of the bulk Lorentz factor and the total 7-ray energy. 
It is, however, interesting to note that Fermi-LAT-detected short GRBs emit higher fluence 
in the LAT range than in the GBM range \20\ [21] which could be due hadronic emission 
processes. 
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4 Ultrahigh-energy cosmic rays from GRBs 



GRBs have long been considered as candidate sources to accelerate UHECRs |22t • The 
energy of protons with 7^ w ^'sat,pi if they were to escape from the GRB blast wave, is 2 x 
lO2Or3/^(0/lO)5^ eV, so GRB 080916C can in principle accelerate UHECRs. The Larmor 
timescale t'^,^{s) ^ 0.2/ y^{^/W)B'^^ at 7p ~ l'sat,p i^ much less than the light-crossing time 
t[Js) ^ 200t^(s), so that escape depends on transport in the jet plasma magnetic field 
impedes escape [23]. Photohadronic processes can assist escape by converting protons to 
neutrons, but are more important when the internal photon density is high, or when the 
variability timescale is small [25]. When ty{s) < 0.01($gbm/10~^ erg cm^'^s"^) / (V^epk) , 
then photopion losses are important in GRB 080916C (see Ref. f26] for applications to GRB 
090510). The shortest variabihty timescale for GRB 080916C observed with INTEGRAL 
is 0.1 s [2], so unless the corresponding size scale for the radiating region was even shorter, 
we can neglect photohadronic processes. But a shorter variability timescale would require 
a larger bulk Lorentz factor F in order to explain detection of multi-GeV photons, in which 
case photohadronic efficiency is reduced [TO] . 

We estimate the rate of long-duration GRBs as energetic as GRB 0809 16C within the 
~ 100 Mpc clustering radius for UHECRs observed with Auger |27|l28j. For maximum total 
energy releases of lO^"^ erg, the GRB 080916C jet opening angle 9j < 100/F = O.I/F3. 
The inferred GRB rate of ~ 2fh Gpc^^yr^^ ^29j at the typical redshift z =1 - 2 is a fac- 
tor ~ 1/10 smaller at w lOOdioo Mpc due to the star formation rate factor and a factor 
fb > 200 larger due to a beaming factor. A 60-^60 EeV UHECR is deflected by an angle 

-1 x/2 1/2 

« 4.4:° ZEnG^QQ '^loo'^i intergalactic magnetic field with mean strength Bnc nG and 
coherence length of Ai Mpc. Defiection causes dispersion in time of arrival of UHECRs [3QJ 
and increases the apparent rate. The corresponding number of GRB sources within « 100 
Mpc with jets pointing within 4° of our line-of-sight is ^ 30{fb / 200) B^^Eq^ xl^'^ . Com- 
plications arising from nonuniform magnetic field geometries (e.g. Refs. |31l I32j ) can lead 
to different values in the rate estimates, but still allow GRBs as UHECR hosts. Thus if 
typical long duration GRBs have a narrow, highly relativistic core accelerating UHECRs, 
then long-duration GRBs could account for the Auger events within the GZK radius. 

5 Discussion 

We have developed a hadronic model based on synchrotron radiation by protons in a highly 
magnetized shocked plasma, to explain delayed onset of high energy (> 100 MeV) emission 
observed with Fermi from GRB 0809 16C. This proton synchrotron spectral component, in 
addition to the electron synchrotron and photospheric emission [33] that dominates in the 
keV - MeV range, initially starts at much higher energy and later sweeps into the Fermi 
LAT range, thus causing a time delay in the prompt phase. Our model is compatible with 
an internal shocks scenario and could also be consistent with an origin of the high-energy 
emission in GRB 080916C from an external shock The late time extended emission 
observed in GRB 0809 16C and other GRBs could be due to the slower cooling of the protons 
in the early afterglow related to the external shock [6l [9] . 

Our approach differs from the conclusions of [35] that the jet outfiow energy in GRB 
0809 16C is dominated by Foynting flux rather than particle energy. Their conclusion de- 
pends on the assumption that the engine radius rg ~ cty/{l + z), whereas ro might be much 
smaller. In our scenario, both the relativistic outflow and the shell on which it impacts are 
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particle energy dominated; only the shocked fluid is highly magnetized. The GBM emission 
could be nonthermal electron synchrotron radiation that includes photospheric emission 
made at smaller radii than the nonthermal synchrotron radiation. Compton cooling of 
shocked electrons is suppressed in highly magnetized shocked plasma, but could be impor- 
tant in a leptonic model with external cocoon photons [36]. Our model for GRB 080916C 
also applies to other Fermi LAT GRBs with comparable T factors and small beaming angles. 

Besides GRBs, UHECRs might also be accelerated in other systems of relativistic out- 
flows, including low luminosity GRBs ^33138], radio galaxies |39j and blazars. UHECRs 
could be formed through neutron escape when photopion processes are important, which 
will require IceCube neutrino detections [lO] to establish. In GRB 080916C, where multi- 
GeV radiation is observed with Fermi, we have shown that synchrotron radiation from 
ultra-high energy protons accelerated in GRB jets explains the delay of the > 100 MeV 
LAT emission with respect to the keV - MeV GBM emission, and that long duration GRBs 
are possible sites of UHECR acceleration. 
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